Nonlocal vortex motion in mesoscopic amorphous Nbo.rGeo.a structures 
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We study nonlocal vortex transport in mesoscopic amorphous Nbo.rGeo.s samples. A dc current 
I is passed through a wire connected via a perpendicular channel, of a length L = 2 — 5 /^m, with a 
pair of voltage probes where a nonlocal response V„i oc / is measured. The maximum of _R„; = Vni/I 
for a given temperature occurs at an L- independent magnetic field and is proportional to 1/ L. The 
results are interpreted in terms of the dissipative vortex motion along the channel driven by a remote 
current, and can be understood in terms of a simple model. 

PACS numbers: 74.78.Na, 74.78.Db, 74.25. Qt, 74.25.Fy 



In a pioneering work Giaver measured a magnetic-flux- 
transformer effect in type II superconductorsi He ap- 
plied a magnetic field B perpendicularly to a sample 
comprising two superconducting sheets separated by a 
thin insulator, passed a current / through one of the su- 
perconductors and measured a voltage developed over 
the other one - where no current was flowing. The in- 
duced voltage was a consequence of an electromagnetic 
coupling of vortices in the two layers. In their recent 
experiment Grigorieva et al^ demonstrated a comple- 
mentary flux-transformer phenomenon associated with 
vortices. They produced mesoscopic amorphous MoGe 
structures of a double-cross shape, consisting of two par- 
allel wires connected at a right angle by a channel of 
a width w = 0.07 — 2 /itm and a length L = 0.5 — 12 
fj,Ta. In a perpendicular B and with / through one of 
the parallel wires a nonlocal voltage Vni appeared over 
the second, current-free wire. This novel, transversal 
flux-transformer effect originated in the in-plane vortex- 
vortex repulsion which conveyed the driving force from 
the current-carrying wire to the vortices in the channel. 
The effect disappeared not only for L exceeding 6-7 /im 
but also for w larger than ^0.5 — 1 fim. When w was 
sufficiently small the force on the vortices in the chan- 
nel was transferred over many intervortex distances, and, 
moreover, Vni was proportional to /. The efficiency of 
transversal flux-transformer effect can be quantified by a 
nonlocal resistance Rni = Vni /I- 

In the experiment of Grigorieva et al^ the local mixed- 
state dissipation was characterised on separate mm-sized 
films, whereas Vni was measured by a low-frequency ac 
method during B sweeps at constant temperatures T. 
An ac method was used because Vni was in nV range, 
i.e. Rni < 5 mf], thus being too small for a dc detec- 
tion. In our work we focused on dc probing of transversal 
flux-transformer effect and measuring Vni ^^nd the local 
voltage V; on a same sample, which was possible in multi- 
terminal amorphous (a-)Nbo.7Geo.3 structures of the ge- 
ometry shown in the inset to Fig. 1. The weak pinning, 
characteristic of the a-Nbo.yGeo.a material used, resulted 
in a dc-measurable Vni and Rni ~ 1 even at very 



low temperatures. The measured nonlocal resistance was 
hence two orders of magnitude larger than in Ref. 

We measured the transversal flux-transformer effect in 
samples of different length by isothermal sweeps of B 
at different applied / and for 0.15 < T < 0.95 T^, 
where Tc is the superconducting transition temperature. 
Vni depends linearly on / in the range / = 0.1 — 1 ^A. 
With increasing B, Rni{B) first acquires a nonzero value 
at a certain magnetic field Bd, then has a maximum at 
some Bp and gradually vanishes close to the upper critical 
magnetic field Bc2- The main representatives of Rni, i.e. 
Bp and Rp — Rni{Bp), behave differently with respect to 
the channel length. Bp is independent of L whereas Rp cx 
1/L, suggesting a vortex velocity Uni oc 1/L at the non- 
local voltage probes. This we relate to the total frictional 
force on the vortices in the channel being proportional to 
L. 

We investigated two structures of the type shown in the 
inset to Fig. 1, where we also assign numbers to the leads 
and define the coordinate system (with the unit vectors 
X, y and z). As in our previous studies;^ the samples 
were produced by combining electron-beam lithography 
with magnetron sputtering but the film thickness was 
increased from 20 nm to d = 60 nm in order to safely 
avoid inhomogeneities detected in samples of a cross sec- 
tion smaller than ^ 5000 nm^. This way the pinning 
was enhanced, as we inferred from the local dissipation, 
but still remained weak enough to permit a nonlocal vor- 
tex motion over several microns. Contacts and 1 are 
used for applying a dc / = Iqi through the horizontal 
wire, exerting a local force on vortices in the y direc- 
tion (since B — —Bz). Having a velocity u = uy, vor- 
tices induce an electric field E — Ek. Combinations of 
contacts i,j = 6 — 9 are used for measuring the local 
voltage drop V. The local vortex pressure is transferred 
from the horizontal wire along the zero-/ channels (the 
length L is indicated for the left channel) toward the 
crosses contacted by leads 2,3 and 4,5, where Vni is mea- 
sured. For the first sample the channel lengths were 2 
and 3 /im, and for the second 3 and 5 /im, so we covered 
the range L = 2 — 5 fim and had two L = 3 fim samples 
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FIG. 1: Local (lines, left-hand scale) and nonlocal (symbols, 
right-hand scale) response for the L — 2 fim channel at T = 
1.4 K {Bc2 = 2.95 T) and 7 = 7oi = 300,600,800 nA. Inset: 
A photograph of the sample, with the designation of the leads 
and the definition of the coordinate system. 



for consistency check. All wires of a single sample had 
the same width, w = 275 nm for the sample with the 
longer channels and w = 250 nm for the other one. This 
small difference did not affect the results presented hence- 
forth. Since dc measurements are invariably burdened by 
sub-/iV parasitic signals, in order to determine Vni{I) for 
a given measurement with leads i,j we recorded both 
Vni,ij{I) and Vni,ij{I = 0) taken at the same i?-sweep 
rate and direction. Vni,ij{I = 0) was different for dif- 
ferent pairs of leads and depended very weakly on the 
B-sweep rate. The background-free nonlocal voltage was 
extracted as Vni{I) = Vni,ij{I) - Vni,ij{I = 0). By this 
procedure we found a very regular behaviour of Vni{I)^ 
which, in particular, for the two L — Z fim. samples 
agreed within the error bars. From the local voltage re- 
sponse between contacts (8,9) we characterised the sam- 
ples in the same way as in our previous work;^ obtaining 
Tc = 2.95 K, the normal-state resistivity pn = 3.7 /ifim, 
-idBc2/dT)T=T, = 2.17 T/K, and the Ginzburg-Landau 
parameters k — 100, ^(0) = 7.15 nm and A(0) = 1.18 /im. 

Results typical of Vni are shown in Figs. 1 and 2. 
The data in Fig. 1 were obtained for the i = 2 /im 
channel at T = 1.4 K {Bc2 = 2.95 T) and with / = 
loi = 300, 600, 800 nA. The lines (left-hand scale) cor- 
respond to the local response Vsg/Joi, and the symbols 
(right-hand scale) to Vni/I = [^45(101) - ¥45(0)]/ hi ■ 
As expected from previous investigations^ of V{I) of a- 
Nbo.yGeo.a, for the given current density J 20 — 50 
MA/m^ the local response at low T depends on /. Note- 
worthy, in contrast to d = 20 nm samples with a weaker 
pinning2i4 the onset of Vi at a certain B = Bdi is essen- 
tially independent of /. Bdi nearly coincides with B^ for 
L — 2 /im, while B^ for L = 3, 5 /zm is higher for 5 % 
(high T) to 20 % (low T) and mutually indistinguish- 
able. Contrary to Vi{I), Vniil) is for the given range of 




FIG. 2: Rni{B) (left-hand scale) for an L = 3 /im channel at 
T = 0.8, 1.6, 2.6 K and / = 600 nA, with the voltage resolu- 
tion shown by the right-hand scale. The curves are offset for 



clarity, 
arrows. 



The meaning of Bd, Bp and Rp is indicated by the 



/ linear for all T.^ 

In Fig. 2 we show how i?„; = Vni/I changes with tem- 
perature. The results refer to an L = 3 /im channel, 
/ = 600 nA and T = 0.8, 1.6, 2.6 K. As in Fig. 1, i?„i = 
up to -B = Bd after which it displays a relatively broad 
peak around B = Bp, which defines Rp — Rni{Bp). Close 
to the Bc2i Rni drops to zero again. While the nonlocal 
resistance in Ref. 2 decreased with decreasing T, vanish- 
ing at T/Tc 0.6, in our experiment we find a nonmono- 
tonic variation of Rni{T). Rni is finite even at the lowest 
measurement temperature (T — 0.4 K, T/Tc ~ 0.14), 
decreases at low T and increases close to Tc with increas- 
ing T (see Fig. 3(b) later). The Rni{B) traces are nearly 
symmetric around 5 = but, especially at low T and/or 
high I, Rp ior B < and B > Q may differ to some 
extent, as seen, e.g., in Fig. 1 for / = 800 nA (1.4 K) 
and in Fig. 2 for T = 0.8 K (600 nA). These differences 
do not necessarily appear only because of variations in 
measurement conditions during the long B sweeps (e.g., 
slightly different bath temperature), they may also origi- 
nate in the Nernst effect^ due to a heating in the current- 
carrying wire.'' However, the Nernst effect changes sign 
by reversal of B and can be canceled by the averaging 
RniiB) = [Rni{+B) + Rni{-B)]/2. Moreover, for our 
samples Rni{+I) — Rni{~I), which rules out a rectify- 
ing effect proposed recently.^ 

Essential information on the nonlocal vortex motion 
is contained in the characteristic quantities Bp and Rp. 
In Fig. 3(a) we plot Bp{T) for all three channel lengths, 
together with 5^2 (T), Bdi{T) and Bd{T) for the L = 2 
/im channel. Bdi and Bd are relatively low and become 
immeasurably small above 2.6 K (this applies to Bd for 
L = 3, 5 /im as well). It can be seen that Bp is inde- 
pendent of L. On the other hand, Rp does depend on 
L, which is shown in Fig. 3(b). RpL vs T, plotted by 
the symbols, exhibits a reasonably-well-defined scaling 
behaviour implying Rp ccl/L. 

We shall now discuss a simple model that leads to a 
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FIG. 3: Characteristic parameters of R„i for the three channel 
lengths, (a) Bp{T) together with Bc2{T), Bdi{T) and Bd(T) 
for the L = 2 /jm channel. The dashed lines are guides for 
the eye. (b) LRp vs T. The solid lines represent the results of 
the calculations (using rjBS and rjT, as indicated) explained 
in the text. 



reasonable explanation of the presented experimental re- 
sults. This model assumes equilibrium between the driv- 
ing force exerted by / on the vortices in the lower cross 
and the frictional damping force on the vortices in the 
channel, as well as the presence of surface barriers enforc- 
ing a dominating y component of u inside the channel. In 
a first step we neglect the pinning and assume that the 
frictional force is linear in u with a velocity-independent 
friction coefficient 77. The effects of pinning are addressed 
later. 

At a vortex density = B/(j)o, where is the 
magnetic-flux quantum, in total n^uP' vortices in the 
lower cross - each experiencing a force J4>Qd, apply a 
pressure p = riipcjjQl /d on the vortices in the channel. The 
corresponding pushing force (per unit vortex length) pw 
is balanced by the force required to move n^Lw vortices 
along the channel against the frictional damping rjUni 
per vortex. This gives Uni = (pQl/riLd cx 1/L. Using 
Vni = wBuni we find 
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The above result holds if the surface barriers are strong 
enough to confine the vortex motion within the channel. 
As shown experimentally in Ref|2, this is not satisfied 
for large w. In any case, surface barriers are essential 
to preserve the uniaxial character of the nonlocal vor- 
tex motion. An important source of the surface barriers 



are the Meissner currents Jm flowing along the channel 
edges and providing an inward-pointing force Fin oc Jm- 
The surface barriers weaken by approaching Bc2 irrespec- 
tively of their exact origin, which my explain the fact that 
Rni{B) does not increase all the way up to B = Bc2 be- 
fore dropping to zero in the normal state. This restricts 
the range of the applicability of Eq. (1) to i? not too close 
to Bc2- On the other hand, Eq. (1) correctly reproduces^ 
Rp cx 1/L and, as we show below, accounts for Rp{T) 
quantitatively in conditions of insignificant pinning. 

For T > 2.7 K the Vi{I) curves are linear beyond 
any doubt, implying a negligible pinning and, more- 
over, 77 « ^7/ of the viscous drag in pure fiux fiow. 
There are two possible dissipation mechanisms that de- 
termine the flux-flow viscosity 77/.^ One is related to Joule 
heating of normal electrons in vortex cores by the E 
therein, as described by the Bardeen-Stephen model giv- 
ing r]f = rjBS = 4'aBc2/ Pn-^^ The other approach, pro- 
posed by Tinkham, attributes the dissipation to a loss of 
the superconducting Gibbs free-energy density Gs{B,T) 
as vortices move and cause depairing and recombination 
of Cooper pairs Ji This process is affected by the time 
T = (^/Ao)[l + {T/Tc)^]/[1 - {T/Tcf] of establishing a 
superconducting state and results in rjf = rjx = 2tttGs 
(Aq « 1.76kBTc is the superconducting gap at T = 0, 
and h the Planck constant). For B not too close to 
Bc2, i.e. in the regime of the applicability of Eq. (1), Gs 
can be approximated by the superconducting condensa- 
tion energy UsiT) « B2^(T)/4k2^o, where no = 47r/10^ 
H/m.'^ Since we know Bc2{T) we can determine both 7755 
and rjT for our samples, which permits a quantitative 
comparison of the experiment and the model. 

At r > 2.7 K the measured Bp « 60 mT < Bc2 is 
fairly constant, so the temperature dependence of the 
corresponding Rp is dominated by that of 77. By calcu- 
lating rjBS a-nd rjT without any adjustable parameter and 
using Eq. (1) we obtain LRp shown in Fig. 3(b) by the 
solid lines, as indicated. Rp{T) is well reproduced in both 
cases. In particular, the steep increase of Rp near can 
now be traced back to the rapid decay of 77/ close to Tc- 

So far we have not included any effect of pinning into 
our analysis, which is justified only in a narrow T range 
close to Tc- At lower temperatures the use of 77/ and the 
measured values of Bp in Eq. (1) results in a rapid in- 
crease of Rni far above the observed values. In addition, 
nonlinearities in V/(/) are observed below 2.7 K, indicat- 
ing that the pinning is no longer marginal. However, the 
fact that the linearity of Vni{I) and Rp oc 1/L are pre- 
served in this regime brings up a possibility of extending 
our model to lower temperatures. This is achievable if 
even in the presence of pinning p remains proportional 
to / and 77 independent of u. Below we argue that these 
properties are consistent with a plausible picture of the 
distribution of the pinning force in our samples. 

We recall that the onset of Vi at Bdi does not depend on 
/, which suggests a critical vortex density for triggering 
the dissipation. This implies that for B < Bdi the flux 
is first trapped in the most-strong-pinning regions of the 
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sampleiiSi When these are saturated a,t B = Bdi vortices 
enter the lower-pinning regions in between - where they 
can move more easily. Only after this saturation the dis- 
sipation starts by the vortex motion along "vortex rivers" 
between the strong-pinning sites/areas. In the current- 
carrying wire vortices in the "rivers" shear plastically 
with the edges, thus enhancing the driving force above 
J0O in depinning the vortices. This process may result 
in a nonlinear V/(/) such as that observed in Fig. 1 but 
cannot propagate significantly into the channel, where 
the external driving current decays exponentially. The 
pushing force is thus conveyed along the "rivers" set by 
the equilibrium pinning properties, and the fraction of 
the contributing vortices is f ^ {B — Bdi)/B. The vor- 
tex pressure is therefore reduced by the same factor /. 

The frictional damping of the vortex motion in " vor- 
tex rivers" has been investigated in samples with arti- 
ficial easy-flow channels embedded in a strong-pinning 
medium. ^■^ Because of the random pinning landscape in 
our samples no commensurability effects between the 
moving and immobile vortex regions are expected and 
the vortex velocity in the "rivers" should respond linearly 
to the applied force. This may explain the observed lin- 
earity of Vni{I), while the magnitude of Rni is reduced 
by the ratio / that accounts for the number of vortices 



in the weak-pinning regions. Although we lack a man- 
ageable model for calculating -q including the effects of 
pinning, the increase of Rp as T is lowered can be ex- 
plained at least qualitatively as a consequence of Bp{T) 
growing faster than r]{T). 

To conclude, we have investigated a transversal flux- 
transformer effect, which manifests itself in a nonlocal 
flow of vortices in a narrow superconducting channel 
driven by a remote dc current. In our low-pinning a- 
Nbo.yGeo.a the effect appears in more than half of the 
superconducting phase diagram, i.e. everywhere where 
the vortices can be moved easily enough to induce dis- 
sipation at very low currents. The effect is two orders 
of magnitude larger than in previous studies^ and we 
observe a nonmonotonic variation of the maximal nonlo- 
cal resistance with temperature, which can be explained 
by the interplay of vortex density and vortex viscosity. 
Close to Tc the data are in quantitative agreement with 
a simple model which is based on the assumption that the 
vortices behave like a weakly compressible fluid confined 
to the superconducting channel by surface barriers. 
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